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I .  INTRODUCTION 

The  H+O  -*■  OH  + 0 reaction  is  an  endothermic  (~  17  kcal/mole) 

1 L . 
reaction  important  in  combustion.  It  continues  to  appear  in  reaction 

2-4 

schemes  describing  various  combustion  processes. 

Trajectory  calculations  have  been  used  to  calculate  cross  sections, 
specific  rate  constants  and  estimates  of  the  overall  thermal  rate  constant 
for  the  reaction.  A new  potential  surface  (Surface  II)  has  been  used 
which  fits  the  known  properties  of  the  HC^  complex  better  than  two 

previous  surfaces  (Surfaces  I and  l' ) . The  new  surface  is  described  by 
two  adjustable  Sato  parameters  instead  of  one  as  in  the  two  previous 
surfaces. 

Surface  I and  associated  results  were  discussed  in  a previous  paper5 
(called  Paper  I).  Surface  I7  corrects  a problem  which  was  subsequently 
found  in  Surface  I.  Trajectory  results  for  Surface  V are  reported  in 

this  report  along  with  those  for  Surface  II.  The  classical  trajectory 

6 ^ ' 

program  CLASTR  , suitably  modified  for  these  LEPS  surfaces,  was  used  for 

all  the  trajectory  calculations.  Surface  I7  shows  little  difference  in 
the  trajectory  results  relative  to  Surface  I.  It  also  has  similar  potential 
contours  for  all  configurations  compared. 

The  results  of  trajectory  calculations  on  Surface  II  indicate 
enhancement  of  the  reaction  rate  due  to  vibrational  excitation  of  the 
0^,  but  the  magnitude  of  this  effect  is  considerably  less  than  on  the 

previous  single  parameter  surfaces.  This  vibrational  enhancement  of 

7 8 

the  rate  is  in  agreement  with  other  studies  ’ on  endothermic  reactions. 


IT  JANAF  Tables,  edited  by  D.R.  Stull  (Dow  Chemical,  Midland,  MI,  1965- 
1967). 

2.  M.B.  Collset,  D.W.  Naegeli,  and  I.  Glassman,  Sixteenth  Symposium 
on  Combustion  (Combustion  Institute,  Pittsburgh,  1977)  p.  1023. 

3.  T.  Miyauchi,  Y.  Mori,  and  A.  Imamura,  Sixteenth  Symposium  on 
Combustion,  (Combustion  Institute,  Pittsburgh,  1977)  p.1073. 

4.  G.  Dixon-Lewis  and  R.J.  Simpson,  Sixteenth  Symposium  on  Combustion, 
(Combustion  Institute,  Pittsburgh,  1977)  p.  1111. 

5.  A.  Gauss,  Jr.,  J.  Chem.  Phys.  68,  1689  (1978). 

6.  CLASTR  is  program  No.  229  in  the  Quantum  Chemistry  Program  Exchange 
Catalogue  (QCPE,  Indiana  University,  Chemistry  Department). 

7.  D.S.  Perry,  J.C.  Polanyi  and  C.  Woodrow  Wilson,  Jr.,  Chem.  Phys. 

3,  317  (1974). 

8.  R.N.  Porter,  L.B.  Sims,  D.L.  Thompson,  and  L.M.  Raff,  J.  Chem.  Phys. 
58,  2855  (1973). 
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The  results  show  clearly  that  the  most  significant  changes  in  specific 
reaction  rates  between  Surface  II  and  Surfaces  I and  I7  occur  for  the 
lower  vibrational  states  (v=0,l,2)  of  the  0^  whereas  the  higher  states 

(v=3,4,5,6)  show  less  of  an  effect.  Surface  II  is  quite  different  from 
Surfaces  I and  I7  in  the  high  interaction  region  (all  nuclei  close 
together).  That  even  small  differences  in  the  shape  of  the  potential 
in  this  region  can  affect  the  reactivity  in  the  lowest  vibrational  level 

9 10 

has  been  demonstrated  by  other  workers 


1 1 . SURFACE  l7  , DESCRIPTION  AND  TRAJECTORY  RESULTS 

The  trajectory  results  for  Surface  I7  are  summarized  in  Tables  I 
and  II,  similar  to  the  presentation  of  results  in  Paper  I.  As  in  that 
study  trajectories  are  run  in  sets  for  which  the  sum  of  the  relative 
translational  energy  (ED)  plus  the  vibrational  energy  (E  ) is  constant. 

The  majority  of  trajectory  points  were  run  with  rotational  quantum 
number  J=1  (corresponding  to  about  a hundreth  of  a kcal/mole  of  rotation- 
al energy) . A number  of  points  were  run  with  J=21  (about  2 kcal/mole  of 
rotational  energy).  Three  types  of  trajectory  events  are  defined: 
reactive,  nonreactive,  and  complex.  Complex  trajectories  are  long  tra- 
jectories that  the  CLASTR  integrator  cannot  follow;  they  are  not  back 
integrable  (energy  and  angular  momentum  are  conserved) . Both  reactive 
and  non-reactive  events  are  back  integrable.  The  complex  trajectories 
have  been  further  subdivided  (Table  II)  into  reactive  complex  and  non- 
reactive complex  trajectories  according  to  the  decision  of  the  CLASTR 
integrator.  Observe  that  in  addition  to  some  points  for  checking  Surface 
I7  results  against  those  of  Surface  I (Paper  I)  additional  points  were 
run  on  Surface  I7  for  vibrational  levels  v=3,2  and  1.  They  fill  in 
nicely  between  the  v=4  and  v=0  results  of  Surface  I.  (See  Paper  I.) 
Specific  rate  constants  have  been  evaluated  for  the  v=3,2  and  1 levels 
assuming,  as  in  Paper  I,  a constant  average  cross  section  over  all 
translational  energies  down,  of  course,  to  the  classical  threshold.  The 
maximum  impact  parameter  used  for  the  cross  section  calculations  was  3 
A.  These  specific  rate  constants  are  tabulated  in  Table  III  for  two 
temperatures,  1600°K  and  2500°K. 

The  difficulty  with  Surface  I was  that  it  had  an  infinity  in  the 
derivatives  of  the  potential  in  the  linear  configuration  (6=180°). 

Surface  I7  eliminates  this  infinity.  This  infinity  causes  little  dif- 
ference in  the  trajectories  results  between  Surfaces  I and  I7  since  the 
system  very  rarely  approaches  the  linear  or  near  linear  configuration 
even  on  Surface  I7  . Because  of  the  large  potential  hole  the  system 


9.  N.  Sathyamurthy , J.W.  Duff,  C.  Stroud,  and  L.M.  Raff,  J.  Chem.  Phys. 
67,  3563  (1977). 

10.  N.  Sathyamurthy,  R.  Rangarajan,  and  L.M.  Raff,  J.  Chem.  Phys.  64, 
4606  (1976). 
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definitely  prefers  bent  configurations.  Also  Surfaces  I and  l'  are 
nearly  identical  in  well  depth  and  equilibrium  configuration.  The  po- 
tential contours  are  similar  for  all  values  of  the  H-0-0  angle. 
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The  only  difference  in  the  mathematical  form  of  the  two  surfaces  I 
and  l'  is  in  the  Sato  parameter.  The  Sato  parameter  for  Surface  l1  is 


| O.S3(sin2Th)2  - 0.  25 (cosTh) 2+0. 2 (y) 2 


TABLE  I.  Summary  of  the  reactive  and  complex  trajectories  for  Surface/ . 

The  total  number  of  trajectories  for  each  initial  0?  vibra- 
tional level  in  each  energy  set  (ED+E  ,)  are  given  also. 

K V 


Ed+E  0- 

R v 2 

kcal  Vibrational 

mole  . , 

Level 

Reactions 

Complexes 

Total 

Traj  ectories 

33 

6 

J=1  RESULTS 

16 

16 

300 

3 

1 

6 

500 

29 

2 

1 

0 

500 

1 

0 

0 

500 

24  { 

4 

7 

24 

300 

. 3 

4 

15 

500 

22  ] 

f 2 

1 

10 

500 

\ 1 

0 

1 

500 

3 

1 

72 

500 

18 

2 

0 

18 

500 

1 

1 

3 

500 

14 

f 2 

0 

44 

500 

U 

0 

18 

500 

33 

6 

J=21  RESULTS 

15 

12 

300 
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TABLE  II. 
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The  complex  trajectories  tabulated  in  Table  I are  divided 
into  reactive  and  non  reactive  groups  using  the  decision 
of  the  CLASTR  integrator. 


Vibrational 

Level 


Total  F 

Complexes C 

J=1  RESULTS 


Reactive 

Complexes 


Non- Reactive 
Complexes 


J=21  RESULTS 


TABLE  III. 


Specific  rate  constants  (at  1600°K  and  2500°K)  for  Surface 
I7  are  given  for  02  initially  in  the  vibrational  levels 
v=3,2,l.  The  populations  of  the  vibrational  levels  relative 
to  the  v=0  level  are  shown.  All  calculations  were  done  with 
an  initial  02  rotational  quantum  number  (J)  of  one.  The 
errors  quoted  are  the  Monte  Carlo  statistical  errors. 


Temperature  (°K) 


Vibrational 

Level 


1600°K 


2500°K 


Specific 
Rate  Constant 
3 

cm 

molecule  sec 

(6.7  ± 4.4)  X 10" 1 

(1.3  ± 1.3)  X 10"J 

(0.29  ± 0.29)  X 10' 

(8.3  ± 5.5)  X 10" J 

(1.9  ± 1.9)  X 10 

(0.61  ± 0.61)  X 10' 


Relative 

Level 

Population 


b = 0.6  X 10  cm. 


Th  = Arcsin(r./x  sin  0) 

I Si 


and 


= (r, 


+ r 


rir2 


cos  0) 


See  Fig.  1 for  a pictorial  representation  of  the  parameters  x , Th,  0, 

SL 


V r2’ 


and  r^. 


The  total  rate  constant  may  be  evaluated  using  the  specific  rate 
constant  data  from  Surface  l'  (Table  III)  and  Surface  I (Paper  I).  The 
total  rate  constant  is  given  by  (only  odd  J appear  in  the  sum  for  O2 ) 


where 


K(T)  = l Fbc(v,J)  k 


v.J 


-E 


f j.  ( 2 J + 1 ) e 


v,J/ 


BC 


\Jv 


vJ 

kT 


BC 


is  the  rotational-vibrational  distribution  function  (Boltzmann 


distribution)  and  QJv  is  the  rotational-vibrational  partition  function. 

The  total  rate  constant  evaluated  from  the  trajectory  results  for 
Surfaces  I and  at  1600°K  is  given  by 


Kt(1600°K)  = (3.0  ± 2.0)  X 10 


-13 


3 

cm 


molecule  sec 


This  value  has  been  evaluated  neglecting  the  possible  contribution  of 
reactive  complexes.  The  error  quoted  is  the  Monte  Carlo  statistical 
error.  The  contribution  from  the  v=0  level  has  been  taken  as  zero.  The 

experimental  value  of  Schott11  is 


Ke(1600°K)  = 1.35  X 10 


3 

cm 


molecule  sec. 


At  2500°K  the  rate  constant  from  the  trajectory  results  of  Surfaces  I 
and  V is 


Kt(2500°K)  = (1.7  ± 0.9)  X 10 


-12 


3 

cm 


molecule  sec. 


11.  G.L.  Schott,  Combustion  and  Flame,  21_,  357  (1973). 
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The  experimental  result  of  Schott  is 


Kc (2500°K)  = 5.9  X 10 


-12  cnT 


molecule  sec. 


As  with  the  1600°K  result  the  trajectory  result  is  a minimum  rate  constant 
(neglecting  reactive  complexes).  If  the  decision  of  the  CLASTR  integrator 
in  regard  to  the  reactivity  of  the  complexes  is  assumed  correct, then  both 
trajectory  results  would  be  doubled,  bringing  the  trajectory  results 
within  a factor  of  two  of  experiment. 


III.  SURFACE  II 

Surface  II  contains  two  adjustable  Sato  parameters.  The  general 

12 

LEPS  potential  for  a three  atom  system  has  the  following  form 

V ■ <V  * 0/  * %'  - <»f  2 2 * “3  2 - “i  “2 

- «lS'  - * D2e 


where 


Q; 

D.6 

(/ 

II 

**- 

• rH 

O' 

i 

l+A . 

1 , 
4 (l+A • ) 

(j*4i 

l 

l 

and 

a . 

D.6 

a. 7 = 
1 

l 

l+A . 

l 

4 (l+A. ) 

l 

l 

' v ' 

i = 

1,2,3. 

The  Morse  parameters  for  the  two  diatomics  OH  and  07  are  given  in 
Paper  I. 


In  the  present  case  the  two  Sato  parameters  are:  1)  A=A^=A, 

(OH  parameter)  and  2)  (0^  parameter) . The  Sato  parameter  A is  closely 

similar  to  the  a'  defined  for  Surface  l'  in  the  previous  section.  The 
Sato  parameter  A is  given  by 

/ / v 2 1 / l 


A = 0 


0. 53(sin2Th)2  - 0.25(cosTh)2  + 0.2(tt/2) 


(see  Fig.  1) . ' ' 

12.  J.T.  Muckerman,  J.  Chem.  Phys.  56,  2997  (1972). 


r 


The  other  Sato  parameter  is  given  by 

A2  = | sin(c  • r2)  + ,15  • £sin(c  • r2)  J | 
l ? ? n 2‘) 

* | 0.53(sin2Th)  - 0.2S(cosTh)  + 0.2  (j)  j 


where 


c = W/5.2  X 10"8  cm. 


The  other  parameters,  Th  and  r.,  have  been  defined  previously  (Fig.  1). 

Surface  II  yields  an  equilibrium  position  for  H0-,  in  better  agreement 

with  recent  experimental  results  than  either  Surface  I or  l'  . For 
Surface  II  the  minimum  energy  is  at  r =0.99,  r_=1.30,  0=105°  versus 

^ ^ 13 

experimental  values  of  r^O.985,  r7=1.36,  0=106°.  Surface  II  is  also 

a better  fit  to  the  shape  of  the  potential  minimum  than  either  I or  V . 

The  force  constants  reported  from  the  analysis  of  experimental  results 

are  fQ^  =6.5  mdyne/A,  fQ0=5.8  mdyne/A  and  fQ=  1.088  mdyne/A.14  The 

values  from  the  new  two  parameter  LEPS  Surface  II  are  f_  = 6.1  mdyne/A, 

UH 

friA=5.1  mdyne/A,  and  f =0.65  mdyne/A.  The  new  two  parameter  surface,  like 
uu  u 

the  previous  surfaces,  fits  the  angular  variation  of  the  potential  as  given 

by  Gole  and  Hayes  (Fig.  2)^.  For  9 < 95°  there  is  some  divergence 
between  Surface  II  and  the  Cl  results;  the  LEPS  fitting  function  is  not 
sufficiently  flexible  to  closely  follow  the  ab  initio  data  in  this 
region.  However,  the  Cl  calculation  itself  is  very  limited,  using  only 
minimal  basis  set  SCF  generated  orbitals.  The  differences  between 
Surface  II  and  Surface  I can  be  seen  by  comparing  the  potential  plots 
for  Surface  II  at  0=80,108,140,  and  180°  (Figs.  3,  4,  5,  6)  with  the 
corresponding  plots  in  Paper  I. 


IV.  TRAJECTORY  RESULTS  FOR  SURFACE  II. 

The  trajectory  results  for  Surface  II  are  summarized  in  Tables  IV 
and  V.  The  mode  of  presentation  is  similar  to  that  in  Paper  I and  to  that 
for  Surface  l'  in  Section  II  of  this  article.  Trajectories  were  run  for 

13.  J.F.  Ogilvie,  Canadian  Journal  of  Spectroscopy,  lj),  171  (1973). 

14.  S.  Farantos,  E.C.  Leisegang,  J.N.  Murrell,  K.  Sorbie,  J.J.C. 
Texeira-Dias , and  A.J.C.  Verandas,  Molecular  Physics,  34,  947 
(1977). 


-O  MINIMAL  BASIS  SET 
-A  C.  I.  WAVEFUNCTION 
— SURFACE  I 
--  SURFACE  H 


100°  120°  140°  160°  180° 

ANGLE  ( 9 ) 


Figure  2.  Comparison  (at  r^O^A,  r2=1.23A)  of  the  angular  dependence 

of  the  LEPS  potential  energy  Surface  II  with  the  ab  initio 
surfaces  calculated  by  Gole  and  Hayes. 
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Figure  3.  Contour  diagram  (in  eV)  of  the  LEPS  potential  Surface  II  for 
0=80°.  Contours  are  labelled  in  eV  with  0.5  eV.  spacing. 
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Figure 


5.  Contour  diagram  (in  eV)  of  the  LEPS  potential  Surface  II  for 
0=140  . 
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TABLE  IV.  Displayed  here  are  the  reactive  and  complex  trajectories  for 
Surface  II.  The  total  number  of  trajectories  is  also  given. 


ER+Ev 

°2 

kcal 

Vibrational 

Total 

mole 

Level 

Reactions 

Complexes 

Trajectories 

J=1 

Results 

38 

6 

22 

0 

500 

33 

6 

26 

5 

500 

( 6 

22 

12 

500 

29 

{ 3 

17 

2 

750 

pi 

I 0 

10 

5 

1350 

24 

H 

/ 3 

14 

10 

750 

\ 0 

16 

S 

1350 

H 

18 

1 

i 3 

3 

26 

750 

\ o 

4 

30 

1350 

H 

J = 

21 

Results 

33 

6 

34 

5 

500 

18  . 

i3 

8 

30 

750 

\ 0 

7 

24 

1350 

TABLE  V.  The  complex  trajectories  for  Surface  II  shown  in  Table  IV  are 


here  divided  into  reactive  and  non-reactive  groups  according 
to  the  decision  of  the  CLASTR  integrator. 


three  different  values  of  E +E  (ED,  relative  translational,  + E , vibra- 

K V K V 

tional  energy)  for  each  vibrational  level  v=0,  3 and  6.  Trajectories 

were  run  for  one  value  of  E.  + E for  the  v=4  level.  The  minimum  cross 

K V 

section  data  of  Fig.  7 are  fitted  with  the  following  function 


The  constants  in  this  expression  for  the  various  vibrational  levels  are 
shown  in  Table  VI.  Interpolation  was  used  to  determine  the  constants 
for  the  levels  for  which  there  is  no  data  (v=l,2,5)  or  only  one  datum 
(v=4) . Note  that  the  v=4  trajectory  point  fits  the  interpolated  curve 
well.  The  thresholds  for  each  cross  section  function  of  the  three 
lowest  levels  (v=0,l,2)  are  a few  tenths  of  a kcal/mole  below  the  clas- 
sical thresholds.  From  these  analytical  cross  section  curves  the  specific 
rate  constants  may  be  derived  using  the  usual  expression 

00 

k ( _8kT_  l 2 f s f(T,E  ) dE 
VJ  ( ^A.BC  f Jo  rvJ  R R 

where  v is  the  vibrational  quantum  number  of  02,  J is  the  rotational 

quantum  number  of  0^  (always  odd  either  1 or  21) , T is  the  absolute 

temperature,  is  the  reduced  mass  of  H.C^,  ER  is  the  collision 

energy  (relative  translational  energy),  and  S is  the  cross  section 

rvJ 

given  by  the  function  previously  described.  Also 


f(T,ER) 


The  specific  rate  constants  are  shown  in  Table  VII  at  two  temperatures, 
1600°K  and  2S00°K.  If  rotational  effects  can  be  neglected  and  assuming 
a Boltzmann  distribution  of  vibrational  energy, the  overall  thermal 
rate  constant  can  be  estimated  (assuming  also  vibrational  levels  above 
v=6  make  negligible  contributions). 

In  the  case  of  the  v=6  and  v=0  states  the  J=21  data  do  not  show 
statistically  greater  reactivity  than  the  J=1  data.  For  the  v=6, 

J=1  state  S =(1.47  ± 0.28)A2,  for  v=6,  J=21,  S = (1.92  ±0.32)A2, 

r6, 1 2r6,21 

for  the  v=0,  J=1  state  S = (.0838  ± 0.0418)A  , and  for  the  v=0,  J=21 

r0,l  2 

state  we  have  S.  = (.1466  ± 0.0552)A  . For  the  v=3  state  the  J=21 
r0 , 21 


i 
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TABLE  VI.  The  parameters  for  the  cross  section  fitting  function 

S = S°  + — tan-1  X (E-E  ) are  given  here  for  each 
rvJ  rvJ  < 

vibrational  level.  Parameters  were  adjusted  for  cross 
sections  fits  to  the  v=0,3  and  6 levels  data  only,  the 
parameters  for  other  levels  were  interpolated.  Rotational 
quantum  number  J is  equal  to  one  (J=l) . 


°2 

Vibrational 

Level 

s° 

rvJ 

(J=l) 

(A2) 

(A’2) 

E 

0 

kcal 

mole 

6 

.60 

122.5 

-8.0 

5 

.49 

150 

-3.7 

4 

.388 

189 

0.7 

3 

.30 

245 

5.0 

2 

.228 

322 

9.4 

1 

.176 

418 

1.37 

0 

.15 

49P 

18.0 

X = 1/2  FOR  ALL  LEVELS 


data  show  statistically  somewhat  greater  reactivity  than  the  J=1  data. 

The  results  are  for  the  J=1  state  S =(1.131  ± 0.652)A2  and  for  J=21, 

2 r3,l 

S = (3.016  ± 1.061)A  . In  all  cases  above  the  J=21  data  are  larger 
r3,21 

than  the  J=1  results.  The  total  thermal  rate  constant  at  1600°K  and 
2500°K  will  be  calculated  assuming  no  rotational  dependence  but  the  above 
data  indicate  that  these  total  thermal  rates  may  be  50  to  100%  low  due 
to  rotation.  This  rotational  dependence  should  be  examined  with  more 
care  in  the  future. 

The  total  rate  constant  at  1600° K from  the  trajectory  work  on 
Surface  II  is 

M1600°K>  - (1.45  ±0.49)  X 10'2  mol“ule  sec. 

This  result  is  quite  close  to  the  experimental  value  given  earlier  (see 
Section  II).  It  must  be  kept  in  mind  that  reactive  complexes  have  been 
neglected  in  this  calculation  and  these  will  augment  the  result,  as  most 
likely  will  the  rotational  contribution. 


TABLE  VII  - Specific  rate  constants  (at  1600°K  and  2500°K)  for  Surface  II 
are  given.  These  are  derived  from  the  cross  section  curves 
shown  in  Fig.  7.  Errors  are  Monte  Carlo  errors.  All  cal- 
culations were  done  with  initial  0,  rotational  quantum 

number  J = 1.  Vibrational  level  populations  relative  to  the 
v=0  level  are  also  shown. 


Temperature (°K) 

Vibrational  Specific 

Level  Rate  Constant 

3 

cm 

molecule  sec. 

Relative 

Level 

Population 

6 

(7. 4*1. 5) 

X 

10' 

■ 11 

3.12 

X 

10' 

4 

5 

(5.9±1.5) 

X 

10' 

-11 

1.13 

X 

10' 

3 

4 

(4. 1*1.2) 

X 

10' 

-11 

4.22 

X 

10' 

3 

1600°K 

3 

(2.0±0.7) 

X 
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■ 11 

1.6 

X 

10' 

2 

2 

(6.4±2.2) 

X 

10' 

•12 

6.3 

X 

10' 

2 

1 

(1.8*0. 7) 

X 

10' 

•12 

2.46 

X 

10' 

1 

0 

(5. 1*  1 . 8) 

X 

10' 

•13 

1 

6 

(9. 4*  1 . 9) 

X 

10' 

-11 

5.  71 

X 

10' 

3 
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•11 

1.3 

X 
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2 

2500°K 
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1 

At  2500°K  the  total  rate  constant  is  given  by 

3 


Kt(2500°K)  = (7.9  * 2.4)  X 


10 


-12  cnT 


molecule  sec. 


This  result  is  some  30%  higher  than  the  experimental  result  of  Schott 
quoted  earlier  (Section  II).  The  trajectory  result  is  low  because  of 
the  reasons  stated  above  for  K^,(1600°K)  and  because  of  the  neglect  of 

vibrational  levels  above  v=6  which  will  contribute  more  significantly 
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(CROSS  SECTION,  Az) 


Figure  7.  The  minimum  cross  sections  are  shown  as  a function  of  trans- 
lational energy.  Monte  Carlo  statistical  error  bars  are 
included  at  each  trajectory  point.  The  data  are  fitted  with 

the  functions  S = S°  + — tan  1 | A(E-E  ) ( where  the 
rvJ  rvJ  K I 0 I 

constants  are  displayed  in  Table  VI. 


I 


than  at  1600°K.  The  error  result  quoted  with  each  rate  constant  is  the 
Monte  Carlo  statistical  error  inherent  in  the  trajectory  method. 


V.  DISCUSSIONS  AND  CONCLUSIONS 

By  using  two  Sato  parameters  which  are  functions  of  the  coordinates 
an  improved  LEPS  surface  (Surface  II)  for  H07  has  been  constructed.  The 

equilibrium  geometry  of  this  surface  agrees  more  closely  with  experiment 
than  the  geometry  of  either  Surface  I or  l'  . However,  the  greatest  im- 
provement is  in  the  force  constants  for  the  HO complex  which  are  much 

better  than  those  of  either  Surface  I or  l'  and  match  the  experimental 
values  quite  well. 

As  for  Surfaces  I and  l'  the  trajectory  results  for  Surface  II  show 
significant  enhancement  in  the  H+07  reaction  rate  due  to  vibrational 

excitation  of  the  0^  molecule.  An  enhancement  factor  of  ~ 150  at  1600°K 

and  — 30  at  2500°K  (between  the  v=0  and  v=6  levels)  is  achieved  for 
Surface  II  as  indicated  by  the  specific  rate  constants.  These  factors 
are  not  nearly  so  large  as  those  for  Surfaces  I and  l'  . The  reason  they 
are  not  so  large  is  clear  from  the  specific  reaction  rate  tables.  At 
the  v=6  level  the  specific  rate  constants  (at  1600°  and  2500°)  for 
Surface  II  are  very  nearly  equal  to  those  for  Surface  I.  However,  as 
the  vibrational  level  is  decreased  the  Surface  II  specific  rates  do  not 

9 10 

drop  nearly  so  rapidly.  Other  workers  ’ have  demonstrated  that  small 
changes  in  the  inner  reflective  wall  of  potential  surfaces  can  dramatically 
affect  reactivity  in  the  lowest  vibrational  level.  The  same  sort  of 
vibrational  effect  is  noted  here  but  the  effect  extends  to  levels 
above  the  lowest  level.  However,  tl  ere  are  rather  significant  differences 
in  the  potential  contours  between  Surface  II  and  Surfaces  I and  I ' in  the 
high  interaction  region  in  our  case. 

An  added  benefit  of  the  new  potential  surface  (Surface  II)  is  the 
dramatic  decrease  in  complexes  formed  for  07  initially  in  the  upper 

vibrational  levels  (v=3  and  above).  Thus  the  minimum  cross  sections  and 
specific  rate  constants  quoted  for  these  levels  are  less  uncertain  due 
to  complexes  than  those  rates  in  Paper  I.  Even  if  all  the  complexes 
were  assumed  to  be  reactive  (a  worst  possible  case)  then  the  specific 
rate  constant  for  the  v=3  level  would  be  some  two  (at  2500°K)  or  three 
(at  1600°K)  times  larger  than  the  minimum  rates.  If  a more  reasonable 
assumption  is  made, such  as  taking  the  decision  of  the  CLASTR  integrator 
as  correct  for  the  complex  trajectories, then  the  minimum  rate  constant 
Ls  only  some  10%  low  (at  both  1600°K  and  2500°K) . The  v=6  level  is 
considerably  less  uncertain  than  the  v=3  level.  Even  for  the  worst  case 
(all  complexes  reactive)  for  v=6,  the  minimum  rate  would  be  low  by  some 
25%  at  2500°K  and  10%  at  1600°K.  For  the  v=0  level  there  is  a significant 


increase  in  complexes  for  Surface  II  over  Surface  I but  there  is  also  a 
considerable  increase  in  reactivity.  In  the  v=0  level  if  all  the  complexes 
were  reactive  the  specific  rate  constant  would  be  some  three  (at  2500°K) 
or  four  (at  1600°K)  times  larger  than  the  minimum  rate.  Assuming  here 
that  the  CLASTR  integrator  decision  is  correct  for  the  complex  trajectories, 
then  the  v=0  specific  rate  constant  is  only  some  20  percent  low  at 
2500°K  and  30%  low  at  1600°K.  In  conclusion,  it  can  be  said  in  view  of 
the  above  results  for  the  specific  rate  constants  that  the  total  minimum 
constant  rate  would  be  in  the  worst  case  a factor  of  three  too  low  due 
to  the  complexes.  It  seems  much  more  likely  to  be  low  by  some  20  to  30%. 

It  would  be  nice  to  eliminate  the  complex  problem  altogether  with  an 
improved  integrator  for  the  trajectories.  Since  the  last  paper  (Paper  I) 

15  17 

several  variable  step  integrators  ' have  been  tried  with  no  improvement 
over  the  CLASTR  integrator. 

As  noted  above  the  reactivity  of  at  least  some  of  the  complexes 
causes  the  total  rate  constants  calculated  in  Section  III  to  be  low. 

Also  these  rates  will  probably  be  low  because  of  the  greater  reactivity 
of  higher  rotational  states.  As  noted  in  Section  III  the  higher 
rotational  levels  (data  for  J=21  level)  are  more  reactive  than  the 
lower  levels  (data  for  J=1  level)  perhaps  by  a factor  of  two.  More 
trajectories  should  be  run  at  various  rotational  levels  to  see  if  this 
tendency  to  greater  reactivity  is  real  or  just  due  to  statistical  error. 

Also  the  contribution  of  higher  vibrational  levels  above  v=6  should 
be  investigated  particularly  at  2S00°K.  However,  including  these  high 
vibrational  levels  (v  > 6)  means  that  some  method  for  including  the 


While  it  appears  that  the  total  rate  constants  calculated  from 
trajectories  for  Surface  II  will  be  higher  than  experiment,  the  rate 
constant  from  the  data  for  Surfaces  I and  If  will  be  low.  Total  rate 
constants  estimated  for  Surfaces  I and  l'  are  some  25%  of  the  experimental 
rates  (at  1600°K  and  2500°K)  if  the  complexes  are  neglected.  If  it  is 
assumed  that  the  CLASTR  integrator  gives  at  least  a statistically 
correct  decision  then  these  rates  will  be  50%  of  experiment. 

As  more  accurate  potential  surface  data  become  available  for  this 
reaction  it  should  be  possible  to  further  refine  the  modified  LEPS 
surface.  Better  dynamics  calculations  should  then  be  possible. 
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